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ABSTRACT 

We use the current sample of ~10,000 zCOSMOS spectra of sources selected with /ab < 22.5 to 
define the density field out to z~l, with much greater resolution in the radial dimension than has 
been possible with either photometric redshifts or weak lensing. We apply new algorithms that we 
have developed (ZADE) to incorporate objects not yet observed spectroscopically by modifying their 
photometric redshift probability distributions using the spectroscopic redshifts of nearby galaxies. 
This strategy allows us to probe a broader range of galaxy environments and reduce the Poisson noise 
in the density field. The reconstructed overdensity field of the 10k zCOSMOS galaxies consists of 
cluster-like patterns surrounded by void-like regions, extending up to z ~ 1. Some of these structures 
are very large, spanning the ~ 50 h~^Mpc transverse direction of the COSMOS field and extending 
up to ~ 0.05 in redshift. We present the three dimensional overdensity maps and compare the 
reconstructed overdensity field to the independently identified virialised groups of galaxies and clusters 
detected in the visible and in X-rays. The distribution of the overdense structures is in general well 
traced by these virialised structures. A comparison of the large scale structures in the zCOSMOS data 
and in the mock catalogues reveals an excellent agreement between the fractions of the volume enclosed 
in structures of all sizes above a given overdensity between the data and the mocks in 0.2 < z < 1. 
Subject headings: galaxies: high-redshift — (cosmology:) large-scale structure of universe — surveys 



1. INTRODUCTION 

Although on scales larger than ~ 100 h~^Mpc the universe is homogeneous and isotropic, on smaller scales structures 
in the universe are o rg:anised hierarchically. T he full view of the cosmic structure was revealed for the first time in 
the first CfA-II slice (|de La pparen t et al.lll986l ). with galaxies distributed in bubble like stru ctures surrounding empty 
regions. Later surveys, particularl y the Sloan D ig:ital Sky Survey (SDSS. fVork et aT]l2000l ) and the two-degree Field 
Galaxy Redshift Survey (2dFGRS, Colless et al. 2001), have confirmed the highly complex hierarchical picture of the 
universe at z ~ 0. Our universe exhibits a range of cellular structures, the cosmic- web (Bond et al. 1996), where 
galaxies define the structures which form patterns of dense compact clusters, elongated filaments and sheet-like walls 
outlining large and almost empty voids. While the network of filaments and clumps sitting in the intersection of 
filaments dominate the visual impression of the cosmic web, most of the volume is in the underdense regions (voids). 

The picture emerging from the surveys of the higher redshift universe is limited by small fields and cosmic variance, 
but the richness and complexity of the cosmic web appears (at least visually) to match that of the local universe out to 
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z ~ 1.5 (e.g. lScQville et al]|2QQ7bl : iMarinoni et al Jl2QQ8h. As redshift increases, the peak in the probabihty distribution 
function of the galaxy overde nsity field shifts towards lower over density values, and the fraction of the volume in the 
underdense regions decreases (|Marinoni eEaD[2008). 

The galaxy distribution is believed to be a (biased) tracer of the underlying smooth and continuous matter density 
field. It has been shown by numerical simulations and by analytical work that the observed large scale structure 
(LSS) pattern is a natural manifestation of the gravitational structure formation process through the amplification 
of density fluctuations and their subsequent collapse. The presence of the cosmic web can be explain ed with the 
tendency of the matter concentrations to co ntract and collapse gravitationally in an anisotropic manner (|Bond et al.l 
[1991 Ivan de Wevgaert fc Bertschingerlll996f ). 

One of the main characteristics of the cosmic-web is its hierarchical nature. The LSS structures cover a broad 
range in (over) density values, physical scales and geometrical shapes. The reconstruction of the density field from an 
observed galaxy sample should ideally preserve all these features. However, the samples of the galaxies are incomplete, 
particularly at high redshifts, and the sampling is often irregular. Different galaxy populations evolve with redshift in 
different ways, and they may trace the underlying mass distribution in different ways at different redshifts. On top of 
this, one has to add an uncertainty in the measured redshifts, magnitudes and other galaxy properties. All of these 
points make the reconstruction of the density field in the Universe non trivial. 

The reconstruction of a continuous, presumably smooth field from a set of measured discrete data points usually 
involves the interpolation and smoothing of the data with some filtering function into a continuous map. Unlike 
the velocity or temperature fields for which we sample the field values in the observed points, for the density field 
reconstruction there are no field samples available. The measured data points (galaxies) are used to define the density 
field itself. A common feature of all the available approaches to reconstruct the density field is that the results of the 
reconstructed density field depend on the method used to carry out the reconstruction. Usually a particular scientific 
application guides the choice of the method to reconstruct the density field, and we briefiy review these here. 

1.1. Applications of the density field 

1.1.1. Extraction of the components of the cosmic web 

An immediate application of the reconstructed density field is the extraction of the four components of the cosmic 
web (clumps, filaments, sheets and voids). While the human eye can easily distinguish between these shapes, there 
does not exist an ideal algorithm yet to extract these features in an algorithmically well-defined way. A lot of progress 
has been made recently, by applying a geometrical classification directly t o the density field, such as the multi-scale 
analysis of the Hessian ma trix of th e density field ( Aragon-Calvo et al.l l2QQ7l ) or the skeleton analysis of the density field 
(iNovikov et al. 2006; Sou sbie et al . 2008) . Alternatively, a classific ation of the co smic web can be done using either a 
linear (|Lee fc Leell2008f ) or non-linear (|Hahn et al.ll2007al )bl: iForero-Romero et alll2008! ) gravitational potential. Some 
of these methods have been applied so far only on numerical simulations and result in the classification of the dark 
matter structures. The reconstruction of the density field from observed galaxy data is based on a very much smaller 
number of objects, and therefore it lacks the spatial resolution of the simulations. However, probably the most difficult 
task to tackle is that the observations provide information on the position of galaxies in redshift space, whi le for th e 
prope r reconstruction of the cosmic web galaxy positions need to be measured in real space (|Lee fc Leell2008l : iLee fc Lil 
I2008D . This is not of course a problem in simulations. 

1.1.2. Correlation between galaxies and environment in which they reside 

Several correlations between galaxy properties and the environments in which they reside are observed in the lo- 
cal Universe. In this context, the environment of a galaxy is given by the number density (or overdensity) of the 
neighbouring galaxies. It has been known for many years that elliptical and lenticular galaxies reside in more dense 
environments than spiral galaxies (e.g. Dressier 1980). Red, older, less star forming galaxies live in mor e dense regions 
than blue, younger galaxies with high star formation rates (e.g. iKauffmann et al.l 120041 : iBlanton et al.l 120061) . These 
trends are observed in environments ranging from clusters to voids. 

Are the early conditions of galaxy formation (the so-called nature scenario) or the environment in which they 
subsequently reside (the so-called nurture scenario), or both of these, responsible for the observed dependences of 
galaxy properties on environment? An obvious way to tackle this question is to go to higher redshifts and to establish 
at which redshift the relations in question are in place. Not long ago, high redshift observations were limited to the 
most dense structures (i.e. clusters) or targeted to detect a specific galaxy population (e.g. Lyman break galaxies or 
Lya emitters). Large and deep redshift surveys are the best probes of the (over-) density field deline ated by g:alaxies at 
all redshifts. Th e recent spectroscopi c high redshift surveys, the DEEP2 Galaxy Redshift Survey (jDavis et al.ll2003D 
and the VVDS (|Le Fevre et al.ll2005f ) and now zCOSMOS (|Lillv et al.ll2007l . 12008 in preph . sample the broad range 
of galaxy population and allow the continuous reconstruction of galaxy environments. Most of the galaxy p roperty- 
envir onment relations observed at z ~ seems to be at place already at z ~ 1 (e.g. Cucciati et al...2006 : Cooper et al.l 
[200I . 

Even in the local Universe, there is an ongoing debate on whether the observed galaxy-environment relations depend 
on the particular s cale a t whi ch the environ ment has been measured. This would potentially give a clue to their origin. 
IKauff mann et al.' (2004) and Bla nton et al.l (|2006l) both find that only the environment measured on the small scales 
(of 1 h~^Mpc) appears to affect the star formation histories of galaxies. In the context of the CDM scenario, this scale 
corresponds to the scale of the individual dark matter halo in which galaxy resides, whereas the larger scale would 
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probe more the location in the cosmic web. 

In the current state of the art of the high redshift surveys, it is very chahenging to reconstruct environments on such 
a smah scale. Galaxies targeted for the observations at higher redshifts are generally more luminous, and their mean 
inter-galaxy separations are thus larger. Furthermore, only a fraction of galaxies is targeted for the observations, and 
this fraction may be substantially smaller than locally. 

1.1.3. Galaxy density field as biased tracer of the matter density field 

In the current cosmological picture it is almost a paradigm that galaxies are a biased tracer of the underlying matter 
distribution. The biasing factor can be a non-trivial function of the scale, of the redshift and of the type of galaxies 
used to reconstruct the density field. Using a statistical approach, one of the ways to infer the biasing function is 
from a comparison of the probability distribution functions (PDF) of the galaxy and matter overdensities (jSigad et al.l 
l2QQQl : [Marinoni et al.| [2QQ5). While the first is based on observations, the PDF of the matter comes from theory and is 
dependent on the assumed cosmological parameters. In the case of the three dimensional overdensities reconstructed 
with a top- hat filter, the PDF of the matter is well described by the log- normal distribution (^Coles fc Jones.. 199L) . 

1.2. Goals of this study 

The main goal of this paper is to reconstruct the galaxy overdensity field in the zCOSMOS region. We reco nstruct the 
overdensity field using the first ~ 10,000 spectra of galaxies from the zCOSMOS survey, up to redshift z ^ 1 (|Lillv et al.l 
12008 in prep ), the so-called 10k sample. High sampling rate and measurements of redshifts with a precision of about 
100 kms"-*^ enable us to delineate the environments of galaxies from the 100 kpc scale of galaxy groups - the scale 
of environment expected to dominate the various galaxy evolutionary processes - up to the 100 Mpc scales of the 
cosmic web. It is one of the major scientific drivers of the zCOSMOS survey to study the role of environment on 
galaxy evolution up to high redshifts. We therefore pay particular atten tion to reconstruct the density field on scales 
that are as small as possible. With this aim, a new method (ZADE, |K ovac et al] lb. in prepi ) has been developed 
to reconstruct the density field using both galaxies with spectroscopic redshifts, and those with only photometric 
redshifts. The photometric probability distribution functions are modified depending on the proximity of galaxies with 
high quality spectroscopic redshifts. We utilise this method to reconstruct the overdensity field at the positions of 
the 10k zCOSMOS galaxies and in any random point in the zCOSMOS volume. The reconstruction has been carried 
out in a number of ways so as to facilitate a wide range of scientific explorations as discussed above. We present the 
three dimensional overdensity maps, and compare the LSS structures in the overdensity field to the independently 
estimated virialis ed structures in the volume of the zCOSMOS survey. Extensive use is made of COSMOS mock 
catalogues (Kitzb ichler fc White! l200"7l ) to estimate the errors on the reconstructed overdensity field and to compare 
the zCOSMOS overdensity field to those obtained from the mock catalogues. Throughout this paper, we assume a 
flat cosmology described with l^rn,o = 0.25 and Hq = 70 km s~^. However, we express the results related to the 
overdensity field using /i, Hq = lOO/i, while stellar masses and absolute magnitudes of galaxies are quoted with an 
explicitly incorporated Hq = 70 km s~^. 

In follow-up pa pers, we study the dependence of star formation properties (|Cucciati et aLllin prepI ) and morphology 
(jTasca et al. in prep) on the environments of the 10k zCOSMOS galaxies presented in this paper^ The envi ron- 
mental d ependence of the luminosity and mass functions is discussed in Zucca et aQ (|in prepI ) and iBolzonella et al.l 
din pred), re s pectively. The biasing function between the zCOSMOS and matter overdensity fields is presented in 
iKovac et al.l (la. in Drepl). O ther st udies cover environmental d ependence of specific gi:alax y populations, such as IR 
galaxies (|Caputi et al.ll2QQ8l ). AGN ([Silverman et aLllsubmittedl ) or post-starburst galaxies (|Vergani et aLllin prepI ). 

2. (GALAXY) DENSITY FIELD RECONSTRUCTION 

A number of methods exist to evaluate the density field, often developed for numerical simulations, e.g. grid based 
methods, smooth particle hydrodynamics (SPH) like methods, Voronoi and Delaunay tessellation field estimators 
(VTFE and DTFE). Here, we concentrate on the problem of reconstruction of the density field from the observational 
data sample of galaxies. 

In practice, reconstruction of the galaxy density field (or galaxy environments) reduces to the (weighted) count of 
objects within some aperture around a set of positions where the density field is to be evaluated. In the general case, 
the density at an observationally defined position r = r(a, ^, z) can be estimated as 

(1) 

The above equation gives the value of a mean density at a position r in the redshift space (averaged over the aperture 
in which the density field is measured), with W a spatial window function, a weight based on the astrophysical 
properties for each galaxy and (j) a function correcting for various observational issues. The difference |r — ri| refers to 
any arbitrarily defined distance between observationally defined points. 
The following points should be noted. 

• The form of equation [T] which we use to define the density at a point is similar to the density estimate in the 
SPH simulations (e.g. iHernquist fc Katzlll98 9), with the exception of the se lection function (/). The density can 
be evaluated at a position of a galaxy (the so-called scatter approach, e.g. IHernquist fc Katd [19891 ) . or at any 
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chosen point in space (the so-cahed gather approach, e.g. iHernquist fc Katzlll989[ ). The points can be of course 
chosen such that they form a regularly spaced grid. 

• The summation in equation [1] goes over those galaxies in the sample that have been chosen to define the density 
field, which we refer to as tracer galaxies. These tracer galaxies might be all galaxies detected in the survey 
(which is commonly fiux limited) or only a subsample of those satisfying some selection criteria, e.g. a criterion 
that is as much as possible independent of redshift, so as to form a "volume limited sample". 

• The function VK(|r — ri|; R) is the kernel used to weight the tracer galaxies, which is a spatial smoothing function 
(e.g. a top-hat or a Gaussian function) and R is the smoothing length. The W function is typically chosen such 
that it weights tracer galaxies depending on their distance |r — ri| from the position where the density is being 
reconstructed. The distance between two points can be defined in various ways. In most of the applications the 
smoothing function W shows some symmetry with respect to the reconstruction point. Commonly, it is of fixed 
shape and normalised to unity, i.e. J WdV = 1. 

• The smoothing length R defines the aperture within which the environment is measured. To obtain a well 
defined value at every point, the field needs to be filtered over a large aperture, such that the shot-noise effects 
are suppressed. However, all properties of the density field on scales smaller than the smoothing filter R will 
be smoothed away. The smoothing length in a given reconstruction procedure can be defined to be of fixed 
or adaptive size. The fixed scale is then the same at every point and should be chosen such that there are 
enough tracer galaxies within the smoothing kernel to be able to reliably reconstruct the density at a given 
point. Unfortunately, in order to keep the number of galaxies in low density regions reasonably large to obtain 
statistically meaningful results, the smoothing length R is then required to be rather large. The statistical 
fiuctuations of the measured density will depend on the number of tracer galaxies within the smoothing aperture, 
i.e. on the density and will be statistically more accurate in the regions with higher density. On the other hand, 
the adaptive scale varies over the volume in which the density field is being reconstructed. It is given by the 
density of the neighbouring tracer galaxies, typically defined by the distance from the point of the density 
reconstruction (which can be a galaxy or any point in space) to the fixed Nth nearest neighbour. The use of 
an adaptive scale leads to a larger dynamical range and higher spatial resolution (in the dense regions) with 
respect to the fixed scale and the purely Poisson noise can also be made constant for a given population of tracer 
galaxies. 

• rrii is the astrophysical weight ("mass") of the tracer galaxy. In the reconstruction of galaxy density fields it 
usually has a value of unity for each tracer galaxies, producing a "number-density". However, rrii can be any 
measured property of the tracer galaxies (e.g. stellar mass). In the ideal case, it would probably be the total 
mass of a tracer galaxy. The resulting "mass" weighted density is supposed to be related to the underlying 
matter density via the bias b factor. However, the bias can be a nonlinear and stochastic function, depending on 
the scale, redshift and type of the tracer galaxies. 

• The reconstructed density field is supposed to be the density field based on the total population of chosen 
tracer galaxies in the observed volume. The function cj) may be introduced to correct for the fact that in reality 
only a fraction of the full population of tracer galaxies is at our disposal when reconstructing the density field. 
The function cj) should account for the observational restrictions such as the non-uniform sampling over the 
observed area, a radial selection function or the redshift success rate of the observed sources in the spectroscopic 
surveys. All of these constraints may depend on the intrinsic properties of galaxies (e.g. the luminosity or the 
morphological type) and as a result it may be very challenging to model all possible dependences into the (j) 
function. 

Unavoidably, the properties of the reconstructed density field will depend on the adopted method. The choice of a 
particular reconstruction method, given by the exact functional forms of W(|r — ri|; i?), (j){T{) and m^, and the choice 
of the tracer galaxies, is guided primarily by the scientific goals, but it will always contain some level of user-specified 
arbitrariness. In particular, one has to optimise between the scale at which the environment is measured and the error 
of the environment reconstructed on that scale, arising because we are using discrete tracers. While it may be desirable 
to measure envir onments on a very small scale, the statistical errors for the smallest scales are also the largest (e.g. 
[Blanton et al. 2003, see also our results on the mock catalogues in Section 4). For some interpretations of the scientific 
results it may be desirable to use a fixed scale to define the environment, because the adaptive scales are smaller in 
the denser regions and larger in empty regions. However, the a daptive scales pro ve to be superior over the fixed scales 
in preserving the complex morphology of the density field (e.g. iPark et ani2007f ). 

One limitation of the outlined scheme is its lack of sensitivity to the geometry of the mass distribution: the used 
kernels are commonly of fixed shape and isotropic. Even smoothing with the adaptive scale will smear out the smallest 
struct ures, particularly the o nes of anisotropic shapes (e.g. Romano-Diaz fc van de Weygaert 2007). Shapiro e t al.l 
(|l996f ) and I Owen et al.l (|l998l ) introduced in simulations an elliptical kernel with axis ratio dependent on the geometrical 
distribution of particles. 

The VTFE and DTFE are other methods to recons t ruct the c ontinuous density fields 
(jBernardeau fc van de Wevgaerd Il996l : ISchaap fc van de Wevgaer^ 120001 : ISchaapl 120071 ). Both VTFE and DTFE 
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are fully adaptive and volume covering methods based on the Voronoi and Delaunay tessellation of the point 
sample, respectively , which divide space into a space-filling network of polyhedral cells (Voronoi tessella tion. iDirichletl 
I185QI : r\^ronoil flQQSl ) or mutually disjunct tetrahedral cells in three dimensions (Delaunay tessellation. lDelonelll934l ) 
according to the local density and geometry of the sampling points. Therefore the VTFE and DTFE methods have 
the advantage over the methods which use kernels of fixed size to better recover the anisotropic structures in density 
fields, such as filaments and walls. 

However, direct applications of the VTFE or DTFE to reconstruct the density field in the current z ^ 1 and higher 
redshift surveys also have some drawbacks. First, the spectroscopic surveys at these redshifts have usually a small 
angular size on the sky (1-2 deg^), and with the current number density of the tracer galaxies the volumes of the 
Voronoi and Delaunay cells will be greatly affected at the edges. Cooper e t al.l (|2QQ5I) used the Voronoi tessellation 
to estimate the local density of galaxies in the DEEP2-like mock catalogues. They conclude that more than 45% 
of the sample galaxies are affected by the edges when using the Voronoi tessellation to reconstruct densities around 
galaxies, compared with ~ 15% when using a cylindrical kernel with R = 1 h~^Mpc. Second, in z ~ 1 spectroscopic 
surveys the fraction of galaxies with a high quality spectrum, which leads to a reliable measure of a redshift used to 
delineate the density field, often is not greater than ~ 50%. It has not been investigated which effects this sampling 
fraction will have on an asymmetric density esti mator such as DTFE. More over, it is not clear how to deal with 
probabilistic objects in the VTFE or DTFE (see Ivan Breukelen eFaH l2QQ6l . who use a Monte Carlo approach to 
sample the photometric redshift probability function of galaxies in the Voronoi tessellation method). The DTFE 
patt erns also contain some artefacts, the most promi nent one being the triangular imprint of the smoothing kernel 
(|Schaapll2QQ7l : [Romang-Diaz & van de Weygaert||2003). 

For some scientific applications, such as biasing, the density field needs to be reconstructed on a fixed scale. In a 
fiux limited survey, the mean separation between galaxies will increase with the redshift for a given smoothing scale 
and increase the shot noise. Even a given population of galaxies will be characterised by some average inter-galaxy 
separation, and therefore with any given galaxy population the density field on scales smaller than this separation will 
be dominated by the shot noise. For fixed apertures, the Wiener filtering technique can be used to deconvolve the 
noise from the reconstructed density field (e.g. Lahav et al. 1994; Hoffman 1994; Zaroubi et al. 1995). This technique is 
based on the minimum variance reconstruction of the density field, which requires the noise model and power spectrum 
of the density field to be known a priori. 

It is common to express the resulting measurement of density as a dimensionless density contrast 6{r) (which we 
interchangeably refer to also as "overdensity" ) defined as 

^(r) = (2) 

where Pm{z) is the mean density at a given redshift. In most of the applications Pm{z) is evaluated as a volume 
average, but it can be estimated also as an average over galaxies (e.g. Cooper et al. 2006). Although we generally 
refer to S also as "overdensity" , we will characterise regions with J > as being overdense, and regions with ^ < as 
being under dense. 

3. ZADE METHODOLOGY 

The photometric redshift technique enables us to obtain redshifts of large number of galaxies in relatively modest 
amount of observing time with respect to spectroscopic surveys. The obtained photometric redshifts can be used 
for a variety of applications, such as luminosity and mass functions. Clearl y, spectroscop i c reds hifts offer a major 
improvement over even a high quality photometric redshifts. For instance, ICooper et al.l (|2005f ) find that for the 
density field reconstruction, if only the photometric redshifts are used, even with uncertainties as small as 0.02(l+z), 
the reconstructed environment in the line-of-sight direction is smeared out on small scales. Only for uncertainties of 
photometric redshifts smaller than 0.005 and when measuring densities in projection, photometric redshift surveys can 
become comparable to spectroscopic surveys in terms of environment reconstruction ([Cooper et al.l 12005). Including 
a properly devised background correction the projected density estimator based on photometric redshifts becomes 
substantially unbiased, even though with a large scatter. This back ground corrected estimator is robust for projected 
densities larger than 10 galaxies per (h~^Mpc)^ (|Guzzo et al.ll2007l ). However, the spectroscopic surveys with redshift 
uncertainties smaller than the velocity dispersion of the small group (i.e. below 200-250 km s~^) are the preferred 
source of data with redshift of sufficient precision to reconstruct the small-scale environments of galaxies. 

To reconstruct the density field on the given scale and with acceptable error one needs to balance between two 
requirements: precise spectroscopic redshifts to accurately trace the overdense and underdense regions, but as large 
number of objects as possible to lower the noise and enable reconstruc tion of the density on the smaller scales. 

We have developed a new algorithm (ZADE, Kovac et al. ' b. in prepi ) which brings together the accuracy of spectro- 
scopic redshifts with the number statistics of the photometric redshifts, in order to broaden the scale of the reliable 
reconstructed environments. Motivated by the literature results of the importance of the scales at which the density is 
measured, our main effort has been made to reliably reconstruct the small scale environments. The algorithm modifies 
the individual redshift probability distributions P(z) that are output by a photometric redshift code (e.g. ZEBRA) 
based on the spectroscopic redshifts of galaxies located nearby on the sky. 

It is well known that galaxies are highly correlated on scales up to 10 h~^Mpc. Therefore, a galaxy that has only a 
photometric redshift probability function P{z)^ is much more likely to lie at some redshifts than at others, depending 
on the accurately determined spectroscopic redshifts of other galaxies that lie near to the same line of sight. 
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As described in more detail in lKovac et alJ (|b, in prepi ), the ZADE algorithm modifies the initial P{z) that is output 
of the photometric redshift programme, by counting the number of objects with spectroscopic redshift within some 
radius (Rzade) as a function of z along the line of sight, to yield a modified probability Pzade{z): 

p _ N{< RZADE,Z)XP{Z) 

^^^^^^"^ - /iV(< RzADE,z)xP{z)dz- 

In Figure ffl we show two examples of the ZADE approach to modify the initial P{z), obtained using the ZEBRA 
code (jPeldmann et aDl2QQ6f ). For the computing purposes, we discretise the ZADE- modified P{z) of every galaxy in 
Az = 0.002 (this does not match the ZEBRA precision). At the current state of the uncertainty in the photometric 
redshifts, there are typically a few structures within the P{z) along a given line of sight (effectively distribution of 
galaxies within Rzade)^ so the ZADE approach is only statistical. 

ZADE has been extensively tested on mock catalogues. We tested our method using five different Rzade values: 
1, 2.5, 5, 7.5 and 10 h~^Mpc, even though in Figure [T] we show results only for three of them (1, 5 and 10 h~^Mpc). 
The value Rzade = 5 h~^Mpc is chosen as an optimal ZADE radius to modify the photometric probability functions 
in o rder to recon s truct the d ensity field without any systematics. 

InE ovac et af] (|b, in prepi) we conclude that the use of the combination of objects with the spectroscopic and ZADE- 
modified photometric redshifts reconstructs the overdensity field at a given scale with smaller errors than the traditional 
approach of using only the spectroscopic sample of objects weighted with some (/)(ri) function in equation [T] (which in 
our tests accounts for a non uniform RA — DEC sampling). Therefore, this enables us to use smaller smoothing scales 
in the reconstruction process, and thereby broadens the dynamical range of the reconstructed environments. 

For every density reconstruction within the zCOSMOS survey we use this combination of galaxies with spectroscopic 
and ZADE modified photometric redshifts. The summation in equation [T] is therefore over all galaxies of a given 
population in the survey volume (and not only over a fraction of galaxies with reliable spectroscopic redshifts) equivalent 
to setting (j){Y\) = 1 for every tracer galaxy. This is one of the major advantages of using galaxies with both spectroscopic 
and photometric redshifts, because in practice it is impossible to model the (/)(ri) function in such a way to take into 
account all the complex selection effects of most spectroscopic surveys. In equation (TJ galaxies with only photometric 
redshift will be counted with their ZADE- modified Pzade{z) according to a smoothing filter W. 

4. RECONSTRUCTING THE DENSITY FIELD WITHIN THE ZCOSMOS VOLUME 

4.1. zCOSMOS survey 

zCOSMOS (Lillv et_alj l2007l . 12008 in p rep) is a redshift survey undertaken in the 1.7 deg^ COSMOS field 



fScoville et al. 2007a). The observations are carried out with the VLT using VIMOS, a multi-slit spectrograph. 
In zCOSMOS-bright a magnitude selection Iab < 22.5 has been applied to select galaxies in the redshift range up to 
^ 1.4. The final sampling rate of the zCOSMOS-bright survey should be a uniform 60-70%. This paper is based on 
the first 10,000 spectra, wh ich have a rather non-uni form sampling pattern, 30% on average [see the Figure with the 
10k zCOSMOS sampling in Lilly et al.' (|2008 in prepI )]. This sample of all spectra yields to the so-called "10k sample" 
of galaxies with only secure redshifts. 

Thanks to the high-quality multi-wavelength photometry of the COSMOS survey (see iTaniguchi et al.l l2007l : 
ICapak et al.ll2007l ). practically all galaxies in the zC OSMOS field can also have a photometric redshift. We employ 
the ZEBRA code in the Maximum-Likelihood mode (|Feldmann et al.ll2006f ) to calculate the full probability distribu- 
tion functions P{z) and photometric redshifts [the maximum likelihood of P{z)] of galaxies in the COSMOS field, 
using typically 10 COSMOS broad bands. The uncertainty in the photometric redshifts that we use in this paper is 
(T2=0.023(l+z) at Iab < 2 2.5. The accuracy in photometric re dshifts will further improve, when new photometric 
data become available fe.g. lllbert et al.ll2008l : ISalvato et al.ll2008i) . 

In our analysis, we complement the "10k sample" of galaxies with ~ 25,000 Iab < 22.5 galaxies which do not have, 
yet, spectroscopic redshifts, but for which photometric redshifts are available. We refer to this composite sample 
as the "10k+30k sample". We additionally refer to the "40k sample" as the set of all galaxies with Iab < 22.5 in 
the area of the zCOS MOS survey. In the analysis below, we will frequently make use of mock catalogues (based on 
iKitzbichler fc Whitell2007 ) for various tests, in which all galaxies in the "40k sample" have accurately known redshifts. 

We use rest-frame absolute magnitudes of galaxies obtained using the ZEBRA code as the best fit template n ormalised 
to each galaxy photometry and best available redshift (spectroscopic or photometric. 'Oesch et al.l [in prepI ). For the 
galaxies with only photometric redshift, a prior which requires Mb to be in the range between -13 and -24 mag is 
imposed. 

Stellar masses are obtained from fitting stellar populatio n synthesis models to the multicolour spectral energy 
distribution (SED) of the observed magn itudes (see [Bo lzonella et al.l [in prepI : [Pozzetti et aLllin prepI ). In this work, 
we use stellar masses calculated usine: the Bruzual fc Ch ariot (2003^" libraries, with the Chabrier initial mass function 
(^Chabrier 2003). The star formation history (SFH) is assumed to decrease exponentially with a time scale r, where 
0.1 < r < 30 Gyr. The ICalzetti et al.l (|200 0) extinction law was used with < Ay < 3 and solar metallicities. The 
final stellar mass is obtained by integrating the SFH over the galaxy age and subtracting from it the so called "return 
fraction" , which is the mass of gas processed by stars and returned to the interstellar medium during their evolution. 
The stellar masses are calculated using the best available redshifts for each galaxy. 

4.2. The choices to measure the density field with zCOSMOS 
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Given the diversity of the scientific apphcations of the density field reviewed above, we have produced a large variety 
of environment measures that are appropriate for different applications. In the following text we discuss the particular 
choices used for the measurement of the zCOSMOS density field. For practical purposes, we use the symbol p for all 
types of densities. 

4.2.1. Tracer galaxies: flux limited or volume limited 

For the 10k zCOSMOS overdensity reconstructions, all galaxies within the zCOSMOS region with Iab < 22.5 are 
used as possible tracer galaxies. Galaxies with high confidence redshift class (fiag > 1.5) and the broad-li ne emitters 
with hi^h con fidence redshift class (fiag > 10+1.5) are counted as one object at the spectroscopic redshift (|Lillv et al.l 
12008 in prepi ). Up to z ~ 1, we do not detect a significant dependence of the assigned redshift class on the type of a 
galaxy. For the rest of galaxies within the zCOSMOS survey limits we use the ZADE-modified ZEBRA photometric 
redshift output. For the few galaxies without a photometric redshift estimate, we use the ZADE-modified N{z) 
distribution of all galaxies in the fields as initial P{z)^ normalised to unity within the redshift interval < z < 1.4 
(only a negligible fraction of galaxies with high confidence redshift class has a redshift z > 1.4). For the few QSOs 
with a low confidence redshift we use the same ZADE-modified N{z) distribution to be their probability function, but 
normalised to fqsOi which is the fraction of the broad-line objects with high-confidence class below z = 1.4 in the 
sample of the broad-line objects with high-confidence class detected at all redshifts. For the 10k zCOSMOS sample, 
fqso = 0.43. 

In summary, we use spectroscopic redshift information of galaxies with a spectroscopic redshift reliability > 99.2% 
(which make about 85% of the spectroscopic sample), combining these with the ZADE-modified photometric redshifts 
for all of the galaxies either without secure redshift or those not yet observed spectroscopically. 

We construct three different samples of tracer galaxies to measure the zCOSMOS over-density field. For the "fiux 
limited sample", all galaxies within the zCOSMOS region with Iab < 22.5 are used as tracer galaxies. In total, there 
are 33,211 such galaxies, of which 8341 have a reliable spectroscopic redshift. In addition, we define two volume 
limited samples of tracer galaxies, which should ideally include the same galaxies at every redshift. Understanding 
of galaxy evolution does not run deep enough to provide a recipe of galaxy evolution which includes all physical 
processes for every single galaxy (e.g. star formation episodes of various activity and merging of galaxies). Therefore, 
to constrain the volume limited samples, we just take into account a passive evolution of galaxies of one magnitude per 
unit redshift interval, which should on average account for most of the evolution for the majority of galaxies. The two 
volume samples satisfy the criteria Mb < —19.3 — z and z < 0.7, and Mb < —20.5 — z and 0.4 < z < 1.0. To include 
galaxies with only photometric redshifts in the volume limited samples we use their maximum likelihood redshift and 
the Mb magnitude at the same redshift. However, in the density field reconstruction, we use the Pzade{z) distribution 
of these galaxies. Galaxies without photometric redshift or Mb are not used to construct the volume limited samples 
of tracer galaxies. The defined volume limited samples are complete samples, and they are selected to be the largest 
possible samples at these redshifts. 

4.2.2. Filter W and smoothing length R 

Starting from the work of lDresslerl (|l980l ). when studying the galaxy properties as a function of their environment, 
it is common to project the tracer galaxies within some velocity interval, and then estimate the surface density. The 
main advantage of using the projected distances instead of the full three dimensional distances is to minimise the 
effects of peculiar velociti es. T he aperture, perpendicular to z, within which the environment is measured can be of 
fixed or adaptive size. Dre sslerl ([l980) for example used the projected distance to the 10th nearest neighbours to define 
a rectangular aperture. 

Following the adopted formalism to describe the density field (equation [Tj) , the filter W in this approach is: 

r^Pi.zADE {z ± 3z) if \rp - r,,p| < D^, \zi - z\<8z {8z = 1000 km s~^) 
otherwise. 

This smoothing kernel is of cylindrical shape, with a radius of Dp in the RA-DEC plane and a length of along the 
line of sight. The distance \rp — rp^i\ is measured between the point where we are reconstructing the density field and 
the individual tracer galaxy projected to the redshift of that point. The Fi^zADE{z ± Sz) is the integrated probability 
of the tracer galaxy to be within the z ^ 5z interval. For galaxies with spectroscopic redshift within this interval 
^i,ZADE{z ± Sz) = 1. We use the interval Sz given by distance of ± 1000 km s~^ to the redshift z of the point where 
the density field is measured. Varying Sz by ± 500 km s~^ around the chosen value of 1000 km s~^ does not make a 
significant difference. 

Using the spatial filter W given in equation [H we measure the density field at the positions of galaxies and on a 
regularly spaced grid. For the former, we include the central galaxy into the counts of objects, if this galaxy passes 
the criterion to be a tracer galaxy. 

We use both fixed and adaptive apertures to estimate the density at a galaxy position. We define the adaptive 
aperture as the projected distance to the Nth nearest neighbour Dn^ excluding the central galaxy from the nearest 
neighbours search, but including it again when estimating the density. For the measurement on a grid we use the 
adaptive projected scale Dp = Dn defined by the distance to the Nth nearest neighbour. Given that we are working 
also with the probability functions (and therefore fractional objects), in practice we measure the adaptive aperture 
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defined by the distance to the Nth nearest neighbour as the minimum distance for which TjiPi^zade^z ± 5z) > N. 
This is valid for both estimates of the density centred on a galaxy and on a grid point. 

Some of the previous studies of the galaxy evolution as a func tion of their environm ent are based on the density field 
reconstructed using the three dimensional Gaussian filter (e.g. ICucciati et aI]l2QQ6l ). To enable a proper comparison 
with such results, we also use a Gaussian filter to reconstruct the density field on scale Rq at the positions of zCOSMOS 
galaxies: 

W{\r-r,\',RG) = I (2^^^)W2)^z,ZAD£;(^z)exp f-^^)'] ^/ \t - t,\ < Rg 
[ otherwise^ 

The |r — Til is a three dimensional separation between positions of a tracer galaxy and a galaxy where the density is 
being evaluated. If the ce ntral galaxy belongs to t he sample of tracer galaxies used, it is also included in the counts. 

For the biasing analysis feovac et aLl la. in prepi ). we use a different scheme to reconstruct the over-density field. We 
want to compare the overdensity field traced by galaxies to the underlying matter overdensity, and therefore we need 
the overdensity field estimated on the same scale as the mass overdensity as available in the literature. Thus we use 
the spherical top-hat filter to reconstruct the density field of galaxies: 

Tx./i I D A / 4 L Pi,ZADE{Zi) if \r-ri\<RTH 

[ otherwise, 

where Rth is the smoothing scale of the top-hat filter. Here, |r — ri| is a three dimensional distance between a tracer 
galaxy and a point where the density is being estimated. An obvious shortcoming of this reconstruction is that it is 
affected by peculiar motions. We include this into consideration for the biasing analysis. 

When the smoothing length of a filter W is defined using the projected distances (as in Equation [3|) we will 
denote the reconstructed overdensity with 6p. In the case when R is defined using the three dimensional distance (as 
in Equations [5] and [6|) , we will denote the reconstructed overdensity with 5. 

4.2.3. Weighting rui of galaxies 

In addition to the number overdensity reconstruction, we also calculate the B-band luminosity and stellar mass M* 
weighted densities. Following equation [TJ these three types of reconstructed environment correspond to the cases of 
rrii = 1, rrii = Lb and = M>^, respectively. We assume that the physical properties of galaxies are constant, i.e. 
they do not change along redshift with PzADEi^i)- Lb and are estimated at the best available redshift for each 
galaxy. 

Weighting each galaxy with its physical property will lead to a galaxy density field still mainly determined by the 
spatial distribution of galaxies, with some additional discrimination based on the properties of galaxies. For example, 
for the same number of galaxies in two identical apertures, the measured number density will obviously be the same, 
while weighting galaxies with their Lb or will add a finer differentiation between the measured densities. 

A shortcoming of the density estimate with 7^ 1 is that an obtained relation between galaxy properties and 
environment can contain some degree of correlation between the studied galaxy property and the property of the galaxy 
used to weight the galaxy counts for the density reconstruction. This is particularly pertinent when including the central 
galaxy itself in the measurement of the environment. For example, the colour-density relation will beco me (artificially) 
more si gnificant when using -weighted overdensities compared to the number overdensities (see ICucciati et al.l 
lin pre pV and the differences betw een the st ellar mass functions in the overdense and underdense regions will appear 
more pronounced (see Bolzonella e t al.lliiTprep ). Nevertheless, rrii ^ 1 weighting can be useful in defining a physically 
better motivated density estimator. 

4.3. Estimating mean density 

To convert the density to an overdensity an estimate of the average density is required. Typically, galaxies in the 
survey itself are used to define the mean density, a value which is supposed to be valid in the whole Universe (at a 
given redshift). Therefore one needs to assume that these surveys are a fair sample of the Universe. 

There are various ways to determine the mean density of a survey. One of the commonly used methods is smoothing 
the observed number distribution of the tracer galaxies along redshift, or fitting this distribution with some assumed 
functional form. Sometimes only the pure number density of objects in some redshift intervals is used, but in this 
case a part of the true variations in the overdensity with redshift will be washed away. The other commonly used 
method is deriving mean densities of galaxies from the galaxy luminosity function, in which case the redshift evolution 
of the luminosity function needs to be modelled. We use here a non-parametric treatment of the luminosity function 
approach. 

We use the full population of tracer galaxies from the zCOSMOS survey (selected from the 10k+30kZADE sample) 
to define the mean volume density as a function of redshift. For each galaxy, we first calculate the maximum volume 
Vmax in which this galaxy could be detected in the zCOSMOS survey, including for completeness an assumed passive 
evolution of Am = Az mag. In the next step the contribution of this galaxy to the density in each redshift bin, up to 
Zmax is calculated as AV{z)/Vmax^ where AV{z) is the volume of the individual redshift bins. Zmax is the maximum 
redshift up to which a galaxy of a given luminosity could be detected in a Iab < 22.5 survey; Vmax is the volume of 
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the zCOSMOS survey up to redshift Zmax- This is carried out for every tracer galaxy and then the contribution of ah 
galaxies is summed to obtain N(z). We present the N(z) distributions for the three samples of the zCOSMOS tracer 
galaxies in Figure [2l Similarly, one can calculate also the luminosity and stellar mass weighted number of objects by 
simply weighting /SV{z)/Vmax contributions by Lb and M*, respectively. 

As desired, the procedure outlined above introduces a smoothing in the very peaky observed distribution of objects, 
as it can be seen in Figure [21 With respect to a more standard smoothing on the observed N(z) distribution, it has 
the advantage of not requiring any assumption on the smoothing length. The kernel in this smoothing procedure is of 
quadrilateral shape, where one side is defined by the distance between Zmin and Zmaxi and the vertical sides are given 
by /Sy{z)/Vmax at Zmin and Zjnax- The remaining side of this kernel is curved, defined by /SV{z) /Vmax in redshifts 
between Zmin and z^ax- The exact shape of the kernel depends on /ab of the individual galaxy. Due to the curved 
side of the kernel, a large number of galaxies at z ~ 0.85 produces a somewhat overestimated N(z) above redshift of 1. 
The produced N(z) distributions using the spectroscopic redshifts or only maximum probability redshifts for the same 
group of objects are almost identical. 

Finally, we obtain the mean density pm{z) by normalising N(z) or its weighted counterparts by the "volume" in 
which N(z) is estimated. When the density is calculated using the three dimensional distances, this is obtained dividing 
N(z) by the volume of the zCOSMOS survey in a redshift bin centred at redshift z of that galaxy (or a grid point). 
To obtain the surface mean density we first integrate N(z) in the interval of ± 1000 km s~^ centred on the redshift of 
a galaxy or a grid point of interest to obtain the corresponding number of objects, which we afterwards divide by the 
zCOSMOS area. 

4.4. Other issues - edge effects 

For all the different types of the density estimators, the effect of the edges of the survey on the estimated density is 
to artificially lower the density, since galaxies are not detected outside of these edges. One of the ways to avoid this 
is to simply exclude all (grid) points at which the aperture size is larger than the distance to the nearest edge of the 
survey. However, in the case of a density estimated using the adaptive aperture, this approach will bias the density 
estimator by excluding regions of the lower density (larger adaptive aperture) in larger proportion than the denser 
regions (smaller adaptive aperture). 

A common way to correct for the edge effects is to scale the measured density with the fraction of the aperture 
which lies within the geometrical limits of the survey, assuming that the tracer galaxies in the part of the cell within 
and outside of the survey limits are distributed in the same manner. 

For the adaptive density estimator the criteria to define the aperture close to the edges of the survey are nevertheless 
still affected by the edge. For example, in the case that the aperture is defined by the distance to the Nth nearest 
neighbour, only galaxies within the survey limits will be used to define this distance. Effectively, for the same density, 
the smoothing length will be larger for the points close to the edges than for the points within the survey. This would 
artificially lower the densities near the edge due to the greater smoothing. We tested an algorithm in which the number 
of required neighbours is lowered close to the edges, according to the fraction of the volume within the survey limits. 
However, based on the tests on the mock catalogues, we find that correcting the densities at the points close to the 
survey limits by simply scaling the density by the fraction of the volume (area) within the survey limits works better, 
and is certainly much easier to implement. 

For all types of the density reconstruction we therefore simply define the edge-corrected density-estimate using the 
following equation: 

Pc = p/f (7) 
where p is estimated using equation [T] and / is the fraction of the adopted aperture that lies within the survey region. 
Depending on the scientific goals, one can also then exclude points with low values of /, if desired. 

5. TESTING THE RECONSTRUCTION METHOD ON MOCK CATALOGUES 

We use mock catalogues to assess the performance of the density field estimator based on the combination of galaxies 
with the spectroscopic and ZADE-modified photometric redshifts when applied on the 10k zCOSMOS survey. The mock 
catalogues for the zCOSMOS survey are extracted from the lightcone-mock catalogues described in Kitzbichler&White 
(2007). These lightcones are based on the dark matter N-body Millenium Simulation run (Springel et al. 2005) with 
ftrn = 0.25, = 0.045, h = 0.73, f^A = 0.75, n = 1 and erg = 0.9. The semi-analytic galaxy modelling that is applied 
to the mock catalogues is described in Kitzbichler&White (2007). Essentially, it follows the model of Croton et al. 
(2006) as updated by de Lucia & Blaizot (2007). The only change is in the dust model, developed to better match 
the observations of the mass functions at higher redshifts (Kitzbichler&White 2007). There are 24 independent mock 
catalogues of an area 1.4° x 1.4° and ^ ^ 7, with a flux limit r < 26 mag. We use 12 of these mocks, randomly chosen, 
to perform various tests. 

The mocks are used in the following way. All mocks are cut to match the exact zCOSMOS area and redshift 
limits. We define the "parent" catalogue (equivalent to the 40k catalogue, Iab < 22.5), in which we assume all tracer 
galaxies in the zCOSMOS volume have a high quality spectroscopic redshift, and the "observed" catalogue, in which 
about 10,000 galaxies have spectroscopic redshifts and the rest of the tracer galaxies have only photometric redshifts 
(analogous to the 10k+30kZADE catalogue in the real data). The distribution of objects with the spectroscopic 
redshifts i n the observed catalogues is modelled (see .Knobel et al...in prep ) to resemble that of the true 10k zCOSMOS 
catalogue (|Lillv et al.ll2008 in prepi ). 
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In the mock catalogues, we model the photometric redshift probability function P(z) as a Gaussian of dispersion 
depending on the selection Iab magnitude. For galaxies with Iab < 22.5 mag we use a = 0.023(l+z) for the purpose 
of the simulations. Galaxies are randomly displaced from their true redshifts by an amount selected from this same 
distribution. We then apply the ZADE algorithm on those galaxies without reliable spectroscopic redshifts in the mock 
catalogues using Rzade = 5 h~^Mpc. We have done a simple check of the robustness of our results with respect to the 
functional form of the probability function by randomly assigning actual probability functions, obtained from ZEBRA 
for the real zCOSMOS galaxies, to the galaxies in the mock catalogue. Our results did not change significantly. The 
shape of the ZADE- modified Pzade{z) is determined mostly by the neighbouring galaxies and not by its initial form. 

For the test on the mock catalogues, the (over) density field is calculated on a grid equally spaced in the angular 
units ARA = ADEC = 2 arcmin and in the redshift Az = 0.002. The density p in the individual grid points is 
estimated using equation [4j with the aperture defined by the distance to the Nth nearest neighbours Dp N- As we 
mention earlier, because we are also work with "fractional" objects, in this reconstruction we count the number of 
objects until SiPi(Az) > N. 

To obtain over densities, the p estimate in the individual grid points is divided by pmeani^)- In this way, the 
reconstructed measure of environment is given in the units of 1 -\- Sp = p/pmean (equation [2]) . For the individual mock 
catalogues, pmeani^) IS obtained from the smooth number distribution N(z) of the full 40k sample of tracer galaxies. 
We obtain the smooth N{z) distributions separately in each mock to closely resemble the uncertainties in pmean{z) 
arising from the limited volume of the real survey. 

We present here only the number over densities, obtained using = 1 for all the tracer galaxies. The tracer galaxies 
are selected from the fiux limited sample of galaxies. In the following plots, the survey area is limited to the central 
0.8x0.8 degrees in right ascension and declination to minimise the edge effects. 

Figure [3] presents a comparison of the 10k-|-30kZADE "observed" and the 40k "parent" overdensities for the indi- 
vidual grid points, using the distance to the 10th nearest neighbour to obtain densities. The green line represents 
the median for the lOkspec+SOkZADE observed, the red lines correspond to the 25th and 75th percentiles for the 
lOkspec+SOkZADE observed overdensities. Binning is done as a function of parent overdensity. Obviously, the ZADE 
method is not perfect - it slightly overestimates the density in the most underdense regions, and it underestimates the 
density in most overdense regions, no doubt because of residual smearing of the initial P{z). However, the reconstructed 
overdensity reproduces the parent overdensity without a significant systematic effect in the majority of regions, within 
the errors. This is shown in Figure [H In all the three panels of Figured! we present a difference between the observed 
10k+30kZADE and the parent 40k overdensities for the individual grid points, plotted as a function of the observed 
overdensity. The density field is reconstructed using three different adaptive apertures defined by the distance to the 
5th, 10th and 20th nearest neighbour. Note that due to the ZADE approach we can compare overdensities obtained 
in the apertures defined by the same N (number of nearest neighbours) in the "parent" and "observed" catalogues 
since they both contain same total number of objects. 

Based on the difference between the 25th and 75th percentiles, the error (statistical scatter) in the reconstructed 
overdensities using the ZADE method is roughly 0.1 — 0.15 dex in 1^6p. It is clear that the error on the reconstructed 
overdensity gets smaller when using a larger number of neighbours, i.e. the density measured at larger apertures. 
Depending on the scientific goal, one has to balance between the smoothing scale and the statistical error on the 
reconstructed overdensity. The presented statistical errors include only those uncertainties arising from the fact that 
we are dealing with a population of galaxies with very different qualities of their measured redshifts (i.e. 10k galaxies 
with spectroscopic redshifts and 30k galaxies with the ZADE modified photometric redshifts). We did not include any 
errors arising from the uncertainties in the spectroscopic redshift or in the measured physical property. Also, these 
errors do not account for any differences between the 40k overdensities and the true matter overdensities. If we redo 
the presented comparison using the full survey area, and not only the inner part, the median and lower and upper 
quartiles in Figures [3] and [4] remains almost unchanged. For example, in Figure [3] there is an indication for the larger 
underestimation of the most overdense regions. 

Figure [3] and Figure [4] are obtained by using one mock catalogue. In Figure [5] we present the difference between the 
observed and parent overdensities as a function of the observed overdensities for all 12 mock catalogues. The errors 
(i.e. the 25th and 75th percentiles of the difference between the observed and parent overdensity) at a given observed 
overdensity are similar for all the individual mock catalogues (and we do not show them in Figure [5]). The errors in 
Figured! can therefore be taken as a reference for the errors in the observed overdensities, when using this particular 
method of the reconstruction of the density field. 

We have chosen to present the errors on the reconstructed density field for this particular method (projected nearest 
neighbour and fiux limited tracer galaxies) as the majority of our studies of galaxy properties are based on the density 
field obtained by using equation [H and the apertures probed by this reconstruction method are smallest. We have 
carried out similar tests to obtain the uncertainties on the reconstructed overdensity when centred on a galaxy, using 
both fixed and adaptive apertures, following equation [H and fiux and volume limited tracer galaxies, weighting them 
with rrii = 1^ rrii = Lb and = . The statistical errors are comparable to the presented errors for the adaptive 
apertures defined by the same number of neighbours as discussed here and the fixed apertures starting from 3 h~^Mpc. 

6. THE lOK ZCOSMOS OVER-DENSITY FIELD 

6.1. Reconstructed environments of the zCOSMOS 10k galaxies 

We present here the reconstructed environments of the 10k zCOSMOS galaxies obtained following the choices 
presented in Subsection 14. 2[ We discuss the obtained dynamical range, the dependence on the chosen aperture and 
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on the weighting function of overdensity values quantified following equation [H where the overdensity is estimated 
centred on a galaxy. The apertures are defined by projecting galaxies within ±1000 km s~^, and we omit direct 
reference to this projection in the following text. 

We obtain the broadest dynamical range of reliably reconstructed local environments of the 10k zCOSMOS sample 
of galaxies by using the apertures defined by the distance to the 5th nearest neighbour in the fiux limited sample of 
tracer galaxies. Based on the tests on the mock catalogue, 5 is the smallest number of neighbours which can be used 
to reliably reconstruct density at all redshifts probed. When using the larger number of objects to define the aperture 
(e.g. 10 or 20), both the most overdense and underdense regions are smoothed out, particularly the most dense ones 
which have a smaller physical size. 

When using the volume limited samples of the tracer galaxies, the distance to the Nth nearest neighbour is equal to 
or larger than the distance when using the fiux limited sample of tracer galaxies. Therefore, the obtained overdensities 
with the volume limited samples of tracer galaxies are also smoothed with respect to the overdensities reconstructed 
with the fiux limited sample of tracer galaxies. 

The distances to the 5th nearest neighbour used to define the aperture for the reconstruction of the density centred 
at a zCOSMOS galaxy for a set of log(l -\-Sp) values are presented in Figure [H for all three samples of tracer galaxies. 
Rescaling this plot one can easily obtain the distances for the other values of since it will vary as {N + 1)^^/^^. 

If we use fixed apertures, the smallest scale at which we are able to reliably reconstruct environments of the 10k 
zCOSMOS galaxies is about 3 h~^Mpc (based on the tests on the mock catalogues), using the density estimate given by 
EquationlH On this and larger scales, the dynamical range of the overdensities is smaller than when using the adaptive 
approach. For example, in the redshift range 0.4 < z < 0.7 where we are complete for both Mb < —19.3 — z and 
Mb < —20.5 — z samples of tracer galaxies, the local environments of the 10k zCOSMOS galaxies are estimated using 
adaptive apertures smaller than or equal to 3 h~^Mpc for 94%, 88% and 62% galaxies for the flux, Mb < —19.3 — z 
and Mb < —20.5 — z volume limited samples of tracer galaxies respectively. The main point to be taken is that when 
using the fixed aperture to measure environments we are not able to differentiate between the most overdense regions 
that can be reconstructed at the adaptive scales. 

The ranges of overdensities discussed above are obtained by weighting tracer galaxies with = 1. There is a good 
correlation between the number overdensities and the Lb and weighted overdensities, because all are primarily 
set by the number of objects. However, the most overdense regions become even more overdense when using both 
rrii = Lb and = M^, allowing even finer differentiation of the most overdense regions. When using the weighted 
counts of galaxies the dynamical range of overdensities is broader than when using the number or Lb weighted galaxy 
counts. As an example, in Figure [71 we compare a subset of differently m^-weighted reconstructed overdensities using 
the Mb < —19.3 — z tracer galaxies in 0.4 < z < 0.7, and aperture defined by the 5th nearest neighbour. Two 
effects are responsible for the observed scatter: the noise of t he mapping between the number of galaxi es and their 
Lb l uminosities or stellar masses M*, and different luminosity (|Zucca et aLllin prepI ) or stellar mass (Bolz onella et al.l 
lin prep) functions in different environments. 

Many of these effects can be seen in Figures [8] and [9l where we present the histogram distributions of overdensities 
reconstructed around the 10k zCOSMOS galaxies. The overdensities are reconstructed in the apertures defined by the 
10th nearest neighbour (Figure [8]) and by the fixed scale of 5 h~^Mpc (Figure [9]), with three types of tracer galaxies: 
Iab < 22.5 (left panels), Mb < —19.3 — z (middle panels) and Mb < —20.5 — z (right panels) and three types of 
weighting: rm = 1 (top row panels), rrii = Lb (middle row panels) and = (bottom row panels). 

The advantage of the use of the volume limited sample is that, at every redshift, we are using the same type of 
objects to define the environment, even though the properties of these tracers can slowly change with redshift due to 
effects of evolution, which we cannot (yet) take out completely. The number density of these tracer galaxies is roughly 
constant, thus the Poisson noise and smoothing scale (if relevant) in the estimated density is on average the same at 
all redshifts. 

On the other hand, when selecting galaxies in the total fiux limited sample, the number of available tracers is larger, 
and one can reliably reconstruct environments on smaller scales than when using the volume limited sample. This can 
be critical for studying galaxy properties as a function of environment. However, different populations of galaxies will 
be used to define the environment at different redshifts and the typical smoothing scale will change with redshift in 
the adaptive approach. If using the fixed scale, systematically a different number of objects will be used to measure 
the environment at different redshifts and therefore the noise component in the estimate will be different at different 
redshifts. Normalisation to the mean overdensity (at a redshift of consideration) can overcome some of these effects. 

6.2. Cosmographical tour of the 10k zCOSMOS survey 

The global picture of the LSS structures traced by the 10k zCOSMOS galaxies is obtained by the reconstruction of 
the density field on the grid filling the zCOSMOS volume. For the zCOSMOS data the overdensity is reconstructed 
on a regular grid with spacing of 1 h~^Mpc in all three directions (along RA, DEC and z axis). All tracer galaxies 
within ± 1000 km s~^ are first set to the redshift of the grid point, and then used to obtain the projected distance to 
the 5th, 10th and 20th nearest neighbour (as for the mocks). The presented overdensity field is the projected surface 
overdensity field Sp^ and the values of the full three dimensional overdensity field S at the same smoothing scale R 
would be larger, 1 + ^ = {1 -\- Sp)^^^'^\ The overdensity field reconstructed in this way is not suited for the cosmological 
analysis, it is more a compilation of local environments as they would be experienced by an arbitrary object residing 
in the zCOSMOS volume. We use this overdensity field for the cosmographical tour of the zCOSMOS survey. 
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We first compare the distributions of overdensities reconstructed with the flux limited tracer galaxies on the grid, 
and at the positions of zCOSMOS galaxies, in Figure [TOl The number of grid points is scaled to match the number 
of zCOSMOS galaxies with a high quality redshift in 0.1 < z < 1. The smoothing of the overdensity structures 
is clearly visible in both the overdense and underdense tails when comparing the overdensity values reconstructed 
using the apertures defined by the 20th with respect to the 5th nearest neighbour, as might be expected. The most 
underdense regions seen on the grid are more underdense than the lowest 1 -\- 6p values obtained centred on galaxies. 
The peak in distribution of overdensities at positions of galaxies is also shifted towards positive overdensity values [by 
log(l + 6^) ~ 0.5 and log(l + 620) ~ 0.4] with respect to the overdensity field reconstructed on the grid points. Again, 
the distribution of overdensit y values on the ^rid d oes not represent the universal volume distribution of overdense 
and underdense regions. In K ovac et al.l (|a, in prep), we will present the distribution of overdensity values obtained 
from the full three dimensional density reconstruction on the fixed scale. 

A visual representation of the zCOSMOS overdensity field is presented in Figures [Til lEl and [131 For these pre- 
sentations we use the overdensity field reconstructed on scales defined by the distance to the projected 5th nearest 
neighbour in the sample of flux limited tracer galaxies. Note that these imply that the actual smoothing scale is 
increasing with redshift (see Figure [6]), and that different populations of galaxies are used in different redshifts to 
reconstruct the overdensity field. However, using the smoothing scale defined by the 5th nearest neighbour in the flux 
limited sample we are able to obtain the broadest possible dyna mical range of the overden sity field and to preserve 
at best the variety of the structures in the zCOSMOS survey. As Hi rauss point out, a flux limited 

redshift survey "...is useful for qualitative and cosmographical description of the structures that are seen, and in some 
sense shows the maximum amount of information in the redshift survey" . 

The reconstructed overdensity field, presented in Figure [Til shows structures in a large range of comoving scales 
covering a spectrum of different overdensities at all redshifts reliably probed by the 10k zCOSMOS survey. Galaxies 
are distributed into cluster-like structures, surrounded by empty, void-like regions up to z = 1. A few points should be 
noted. First, the increased smoothness and extension of the structures in the overdensity field towards higher redshift 
is an artefact of the increased smoothing scale with redshift. Second, for the smallest scale structures the imprint of 
the filter W used is also clearly visible, especially in Figure [TTl Even though sampling in the zCOSMOS field is not 
uniform, the resolution of similar structures is the same, given that we use ZADE to take into account galaxies without 
spectroscopic redshift and keep their position on the sky. 

Figure [El provides another view of the complex structure of this overdensity field, but compressed by a factor of 
~ 3.5 in redshift direction. Here, we present the LSS delineated by the isosurfaces enclosing regions with l^6p > 1.5, 
3, 5 and 10 (1 + (5 > 1.8, 5.2, 11.2 and 31.6) going from the left to the right, respectively. The isosurface structures up 
to 1 + (5p > 5 are connected over the transverse comoving scales covering the full zCOSMOS area, and coherent over 
hundred or more comoving h~^Mpc (or more than Az ~ 0.1) in the radial direction. 

Three large structures dominate the zCOSMOS cosmic web at z < 1. These structures are located at z ~ 0.35 
(comoving distance Dc - 980 h-^Mpc) , - 0.7 {Dc - 1800 h-^Mpc) and - 0.85 {Dc - 2100 h-^Mpc) and they 
correspond to the peaks already visible in the number distribution of galaxies with redshift (see Figure [2]). The 
connectivity of the two largest high redshift structures extend over 200 h~-^Mpc in the radial direction. Even though 
we are working with a flux limited sample of tracer galaxies, meaning that we are detecting only the brightest and 
presumably most massive galaxies, at z ~ 0.9 the smoothing scale for the overdensities 1 -\- Sp > 3 are still below or 
slightly larger than 1 h~^Mpc, see Figure [6l Based on the tests on the mock catalogues, which include the selection 
effects of the 10k zCOSMOS survey, the existence of these huge overdense structures at high z is not an artefact of 
our reconstruction method. The other striking element in the zCOSMOS overdensity field is a very small number of 
1 -\- 5p > 3 (and above) structures in 0.4 ^ z ^ 0.6 (1100 ^ Dc ^ 1580 h~^Mpc) except for an overdense sheet at 
z - 0.53 {Dc - 1420 h-^Mpc). 

Complementary to the large positively overdense structures, we identify in the zCOSMOS 1 -\- Sp field also the 
structures enclosed by the overdensity values of 1 + = 0.15 {1 -\- S = 0.06) and 5p = 0.25 {S = 0.13) which contain 
only regions with 1 -\- 5p lower than these (the underdense regions. Figure [13]). The identified regions are at least 
6.67 and 4 times less dense than the mean density, respectively. These regions with low galactic density do not show 
strong clustering in either the transversal (RA-DEC) or radial (redshift) direction. Also, they appear to be more 
homogeneously distributed over the zCOSMOS redshift range, even though they are less present in the redshift ranges 
of the three most overdense structures, being particularly absent at z ~ 0.35 {Dc ~ 980 h~^Mpc). 

Clearly, the amount of LSS in zCOSMOS field varies with redshift. Quantitatively, this cosmic variance of the 
observed structures is shown in the left panel in Figure [141 where we plot the fraction of the volume of the zCOSMOS 
survey enclosed within the isosurface structures of a given 1 -\- Sp value. We split the sample in four Az = 0.2 slices, 
starting from z = 0.2. The fraction of the volume occupied by the structures with at least 1 -\- Sp values is rapidly 
decreasing with l-\-Sp. The two higher redshift slices are statistically more representative (they occupy larger volumes), 
and we see that for the slices above z = 0.6 about 44% of the volume is in the structures with 1 -\- Sp > 1, and about 
10% of the survey volume is in the structures with Sp > 3. The cosmic variance is particularly noticeable between 
the two lower redshift slices. At 1 -\- Sp = 3 and 1 + (5p = 15, the volume fractions in those structures in the two lower 
redshift slices are different by a factor of ~ 3 and ~ 4, respectively. 

The fraction of galaxies residing in the isosurface overdensity structures is much higher than the volume fraction 
of these structures. Taking the two higher redshift slices as statistically more representative, we conclude that about 
50% of the galaxies reside in structures with 1 + > 3.5, while about 50% of the stellar mass and about 50% of the 
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B-band luminosity at a given z inhabits the structures with 1 + > 4 (three right panels in Figure [T4j). The overdense 
structures are few tens time more important in terms of the fractions of their baryonic content, than in terms of the 
volume fractions which they occupy. 

6.3. Comparison to the other LSS measures 

The high resolution of the zCOSMOS-bright spect ra allows also to identify virialised groups of galaxies with velocity 
dispersion of ~ 250 km s~^. iKnobel et al.l (|in"prepl ) have applied the friends of friends and the Voronoi based group 
finding algorithms to the zCOSMOS galaxies with the high confidence redshift. In total, their optimal group catalogue 
contains 800 groups with at least 2 detected members up to redshift of 1 in the zCOSMOS volume. From these, 151 
groups have at least 4 members detected. Tests on mocks show that the vast majority of all of the identified groups 
will be virialised objects. 

The other tracer of the virialised structures is the hot baryonic gas detected in the X-rays. iFinoguenov et all (|2007h 
carried out the identification of the X-ray clusters in the COSMOS field using the XMM- Newton observations. The 
updated catalogue, improved with new additional X-ray data and usin^ the s pectroscopic zCOSMOS redshifts, contains 
2 18 X-ray clus t ers de tected with high confidence (jFinoguenov et aLllin prepi ). 

iGuzzo et al.l (|2007l ) used projected densities (reconstructed with the photometric redshifts). X-ray surface brightness 
and one of the first weak-lensing convergence maps to describe the extendend structure at z ~ 0.7 in the COSMOS 
field. A comparison of the overdensities reconstructed at the positions of all zCOSMOS galaxies and of those galaxies 
residing in the virialised structures is presented in Figure [151 The majority of galaxies detected in the X-ray clusters 
or the richer optical groups reside in extremely overdense regions. The galaxies defining the poorer groups extend to 
much lower overdensities than the rich > 4 optical groups and X-ray clusters, and they almost completely avoid the 
most overdense regions. 

In addition, in Figure [16] we prese nt the visual com parison of the zCOSMOS overdensity field with the bound 
structures: X -ra y c lusters (left; Finog uenov et al.l lin p rep) and zCOSMOS groups with at least 3 members (right; 
iKnobel et aLllin prepI ). There is an overall good correspondence of the overdense regions and the bound structures in 
the zCOSMOS volume. Almost none virialised structure is detected in regions with 1 + < 3. 

More detailed insight into the spatial distribution of the overdense and virialised structures is presented in the panels 
of Figure \17\ We use the same data as in the previous figure, but now the overdensity field is projected in redshift 
slices of Az = 0.025 width, starting from redshift 0.2. We use the cut of 1 -\- Sp > 6.67 to define the overdense regions, 
higher than in Figure [TH For a contrast, we also plot the equivalently underdense regions defined by isosurfaces with 
l-\-Sp < 0.15. The X-ray structures and the optically defined groups with at least 3 detected members reside inside of 
the LSS defined by the chosen overdensity value in most of the cases in the whole redshift range 0.2 < z < 1. However, 
there are some virialised structures which do not live in these most overdense regions. In fact, as we have already seen 
in Figures [15] and [161 the virialised structures trace also the less overdense regions. Moreover, the apparent overlap 
between some of the virialised structures and underdense regions in Figure [17] is only due to the projection effects. 

6.4. Comparison of the 10k zCOSMOS LSS to the LSS in the mock catalogues 

The overdensity field reconstructed using the galaxies detected in the zCOSMOS survey is highly complex, resembling 
the network of the local cosmic web up to the highest redshift z = 1 probed. We compare here the zCOSMOS 
overdensity field to the overdensity fields in the "10k+30kZADE" mock catalogues. The mock catalogues and the 
density field reconstruction procedure on the mocks were discussed in Section [5] For the exact comparison of the data 
and the mock overdensity field, we reconstruct the overdensity field of the 10k zCOSMOS galaxies on the same grid 
as was used for the mock catalogues, defined by ARA = ADEC = 2 arcmin and Az = 0.002 and using the flux 
limited tracers of galaxies. We did not calculate the edge corrections for all the mock catalogues (which is a very time 
consuming process computationally), and instead compare the overdensity fields uncorrected for edge effects for both 
the data and the mocks. Given that we are only interested in comparing the structures in the overdensity fields in the 
data with structures in the mocks, and not to draw any scientific conclusions from the distributions of the detected 
structures, this should not matter. 

In Figure [iHl we show the three dimensional distribution of the l-\-Sp = 3 isosurface in the overdensity field estimated 
on the angular grid of the 10k zCOSMOS galaxies and of the mock catalogues. The visual inspection of the LSS defined 
by the 1 -\- Sp = 3 isosurfaces leaves the impression that there are more large structures in the real data than in the 
mock catalogues. 

For a quantitative comparison of the overdensity field of the real zCOSMOS and mock catalogues we implement 
a volume filling statistic. We calculate the fraction of the survey volume in which the overdensity value is above a 
given threshold. As a reference, we use the value obtained by averaging the individual statistics from the 12 mock 
catalogues. We use the standard deviation of the 12 mock results as an error estimate, which is dominated by the 
cosmic variance. 

First, we compare the overall distribution of the survey volumes in structures of a given 5p value in the 0.2 < 2: < 1 
redshift range using the overdensity field reconstructed in the apertures defined by the distance to the 10th nearest 
neighbour of the fiux limited tracer galaxies (Figure [191 left). Even though we use the angular grid to reconstruct the 
density field, we count the comoving volumes (in [h~^Mpc]^) of the individual cells to properly calculate these volume 
fractions. There is an excellent agreement, within la errors, between the fractions of the volumes in the isosurface 
structures in the data and in the average mock. 



14 



We carry out the same analysis dividing the redshift range in four intervals: 0.2 < z < 0.4, OA < z < 0.6, 
0.6 < z < 0.8 and 0.8 < z < 1 (Figure \T9\ right). While the overdensity distributions of the data in 0.2 < z < 0.4 and 
0.6 < z < 0.8 redshift ranges are in relatively good agreement with the mock results, the redshift interval 0.4 < z < 0.6 
is underdense in the data with respect to the mocks, while the situation is reversed in the redshift interval 0.8 < z < 1. 
Obviously, at any redshift bin and at a given Sp^ the real data are affected by the cosmic variance, and we can not see 
the growth in the cosmic structure over redshifts, as it is visible in the mock curves. 

We follow up on this with a more detailed comparison of the distribution of "sizes" (volumes) of the structures 
above a given overdensity in the data and in the individual mock catalogues. We calculate the size of a structure 
by adding comoving volumes of the connected grid cells with overdensity above a chosen value. The grid cells need 
to have at least one common side (either in RA, DEC or z direction) to be considered connected. The results are 
presented in Figure [20] for l-\-5p > 3. The individual plots correspond to the survey volume fractions contained within 
the structures of at least the indicated size. It is noticeable that there is not a single mock catalogue which contains 
as much volume as the real data in the largest structures in 0.2 < z < 1. As a check we also recalculate the same 
statistics for the data when using the mean density estimated following the same smoothing scheme as applied to the 
mocks (even though this smoothing is not ideal for the data). In this case, the difference between the data and the 
mocks is even larger. 

Carrying the same analysis in the narrower redshift intervals we realise that the majority of this difference is accounted 
for by the large structure in 0.8 < z < 1. Already, when we limit the redshift range to 0.3 < z < 0.9, the data are not 
so different from the mocks. To conclude, except for the large structure in the highest redshift bin 0.8 < z < 1, the 
fraction of the volume within the structures of a given size in the data and in the mock catalogues is in reasonable 
agreement. At the current state, given that there is a disagreement in only one Az slice, it is difficult to say whether 
the large structure in 0.8 < 2: < 1 reflects only a cosmic variance, or it is an unusual object for the current cosmology. 

7. SUMMARY 

We have used the first ~ 10,000 spectra from the zCOSMOS bright survey to reconstruct the density field in the 
survey volume up to 2; = 1. We use a new method for the reconstruction, which is based on the combination of the 
high quality spectroscopic redshifts and the ZADE-modified photometric redshifts of galaxies without spectroscopic 
redshifts. Our method enables us to reliably reconstruct a broader range of environments than it would be possible 
by using only galaxies with spectroscopic redshifts. 

We use the weighted counts of tracer galaxies within various apertures to measure the galaxy environments in the 
zCOSMOS volume. Given the variety of the scientific applications, we carry out the reconstruction in different ways. 
We construct flux and volume limited samples of tracer galaxies, specify fixed and adaptive apertures (characterised 
by distance to the Nth nearest neighbour), define three dimensional and projected distances. We weight tracer galaxies 
in three different ways: with unity, B-band luminosity and stellar mass. 

We present in detail the density field reconstructed on the grid filling the zCOSMOS volume up to z = 1. The 
apertures are defined by the distance to the Nth nearest neighbour (5th, 10th and 20th) projected within ± 1000 
km s~^ to the redshift of a point where the density is being measured. We use the flux limited tracer galaxies, as 
this sample of tracers allows to reconstruct the density field with the broadest dynamical range and with the greatest 
detail of the structures. The reconstructed zCOSMOS overdensity field consists of cluster-like structures, surrounded 
by void-like regions, showing a complexity of the cosmic web up to 2: ~ 1. 

The regions in the density field enclosed by the 1 + (5p > 3, or higher cutoff values, are well traced by the virialised 
objects in the zCOSMOS volume, i.e. X-ray clusters and optical groups with at least 3 detected members. The galaxies 
defining the poorer groups (with 2 or 3 detected members) are found to live also in much lower overdensities than the 
rich > 4 optical groups and X-ray clusters. 

Further on, we have compared the LSS in the zCOSMOS data and in the mock catalogues. There is an excellent 
agreement of the fractions of the volume enclosed in structures of all sizes above a given overdensity between the 
data and the mocks in 0.2 < 2; < 1. However, for the chosen 1 + = 3 value, there is more volume enclosed in the 
large structures in the data than in any of the used mocks; we want to stress that this difference is mainly driven by 
the existence of a very large structure in the zCOSMOS, centred at 2; ~ 0.9 and extending in radial direction over 
Az - 0.2. 
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0.4 



Fig. 1. — ZADE approach to modify the initial P(z), discretised in Az = 0.002. The photometric redshift probabihty functions are 
presented in the lower panels. The P(z) output from ZEBRA is presented in red and the ZADE-modified PzADsiz) are presented in 
blue (RzADE ^ 1 h~^Mpc), black (Rzade ^ 5 h~-'^Mpc) and magenta(Rz ad E < 10 h~-'^Mpc). The number counts of objects with 
spectroscopic redshifts within Rzade radii at a given z are presented in the upper panels, using the same colour coding as in the lower 
panels. The left panel is for an object with I ab =21.94 mag, the right panel is for an object with Iab=22.38 mag. Note that the number of 
objects with spectroscopic redshift within Rzade < 1 h~-'^Mpc is very small, and therefore N(R < Rzade=^ h~-'^Mpc) is zero at almost 
all z (upper panels). See text for more details. 



120 



100 



80 - 



~1 1 1 T 



~1 1 1 T 



~1 1 T 



number density 




Fig. 2. — Redshift distribution of galaxies in the zCOSMOS area. The continuous lines correspond to the smoothed N(z) distributions 
obtained by weighting galaxies according to their AV(z)/Vmax contribution in Az = 0.002 intervals as described in the text. The dotted 
lines correspond to the histogram distributions of tracer galaxies in the redshift bins of Az = 0.05 and they are divided by 25 to match 
the redshift bin of the smoothed N(z). Magenta represents the flux limited sample of tracer galaxies, blue represents the sample of galaxies 
with Mb < —19.3 — 2; and green represents the sample of tracer galaxies with Mb < —20.5 — z. The last two samples are volume limited 
up to 2; < 0.7 (blue curve) and z < 1 (green curve), respectively. 
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0,1<z<^0.4 0.4<z<=0.7 0.7<z<=1.0 




-2 2-2 2 -2 2 

log(l+<5P,) [N-10] 

Fig. 3. — Comparison between the 10k+30kZADE observed and 40k parent overdensities. The green hne represents the median for the 
lOkspec+SOkZADE observed, the red hnes correspond to the 25th and 75th percentiles for the 10kspec+30kZADE observed. Binning is 
done along the parent overdensity axis. Overdensities are obtained by using the distance to the 10th nearest neighbour. Note that 50% of 
all points are contained between the red lines. The lower indices p and g indicate that the overdensity has been reconstructed using the 
projected distances to define a smoothing length R and has been reconstructed on the grid, respectively. The upper index refers to the 
type of the mock catalogue; "p" stands for the parent and "o" for the observed catalogue. See more details in the text. 
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Fig. 4. — The difference between overdensities reconstructed in the observed catalogue, using the ZADE formalism, and in the parent 
catalogue, plotted as a function of the observed overdensity. The density is obtained by measuring distance to the 5th, 10th and 20th 
nearest neighbour in the top, middle and bottom panel, respectively. The continuous line represents median for the measured difference 
between the two reconstructed overdensities, the dotted lines correspond to the 25th and 75th percentiles for this difference. The indices 
are defined in Figure [3] 
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Fig. 5. — The difference between overdensities reconstructed in the observed catalogue, using the ZADE formahsm, and parent catalogue, 
plotted as a function of the observed overdensity. The density is obtained by measuring distance to the 10th nearest neighbour. Each 
continuous line corresponds to the median of the difference between the two reconstructed overdensities (observed and parent), in the bins of 
the observed overdensity, for a single mock catalogue. Individual panels include results for 12 mock catalogues. The observed overdensities 
are presented in the range — 1.2 < log(l + Sp^g) < 1.2. The indices are defined in Figure [3] 




Fig. 6. — Projected distances, defined by the distance to the 5th nearest neighbour, used to define the aperture for the reconstruction of 
densities centred on the 10k zCOSMOS galaxies. The continuous, long-dashed and short-dashed lines are for the /ab < 22.5, Mb < —19.3—2; 
and Mb < —20.5 — z samples of tracer galaxies. The lines for the overdensities of 0.5, 1, 3 and 7 in 1 + (5p units are presented in green, 
blue, pink and red, respectively. 
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log(l + (55)[MB<-19.3-z] log(l + (55)[MB<-19.3-z] 

Fig. 7. — Comparison of the overdensities measured around the 10k zCOSMOS galaxies using equation |4] with the Mb < —19.3 — z 
volume limited sample of tracer galaxies and within the apertures defined by the distance to the 5th nearest neighbour. We compare the 
B-band luminosity and stellar mass weighted overdensities to the unity weighted overdensities in the left and right panels, respectively. The 
magenta line is the median and the cyan lines are the lower and upper quart iles of the presented distributions, where the binning is carried 
out along the abscissa axis. There is a good correlation between the number overdensities and the Lb and M* weighted overdensities. 



22 




log(l+(5p,,J[I,3<22.5] log(l + (5p,,J[MB<-19.3-z] log( 1 + c5p,, J[Mb<-20.5-z] 

Fig. 8. — Histogram distributions of number of galaxies in the bins of overdensities log(l -\- 6p). The overdensities are reconstructed 
at positions of the 10k zCOSMOS galaxies, using the cylindrical adaptive filter W (equation \^ with smoothing length defined by the 
10th nearest neighbour. The distributions are shown separately for three different redshift bins: z < 0.4 (blue), 0.4 < z < 0.7 (red) and 
0.7 < z < 1 (green). We use mass-weighting rrii = 1, rrii = Lb and rrii = M* for the reconstructions presented in the upper, middle 
and lower panels, respectively. A set of these three distributions of overdensities reconstructed for the /ab < 22.5, Mb < —19.3 — z and 
Mb < —20.5 — z samples of tracer galaxies is presented in the left, middle and right, respectively. 
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Fig. 9. — Histogram distributions of number of galaxies in the bins of overdensities log(l + Sp). The overdensities are reconstructed at 
positions of the 10k zCOSMOS galaxies, using the cylindrical fixed filter W (equation 0} with smoothing length defined to be 5 h~-'^Mpc. 
The distributions are shown separately for three different redshift bins: z < 0.4 (blue), 0.4 < z < 0.7 (red) and 0.7 < z < 1 (green). We use 
mass-weighting rrii = 1, rrii = Lb and rrii = M* for the reconstructions presented in the upper, middle and lower panels, respectively. A 
set of these three distributions of overdensities reconstructed for the /ab < 22.5, Mb < —19.3 — z and Mb < —20.5 — z samples of tracer 
galaxies is presented in the left, middle and right, respectively. 
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Fig. 10. — Histogram of the overdensities in the zCOSMOS volume calculated on the grid points (blue) and at the positions of zCOSMOS 
galaxies (magenta). The histogram of the grid points is scaled to the number of galaxies with high quality redshifts. Galaxies are not 
distributed homogeneously in the zCOSMOS volume, they form clumps of the overdense regions leaving some of the volume empty. 
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Fig. 11. — Cosmographical tour of the zCOSMOS survey. The 10k zCOSMOS overdensity field is reconstructed on the grid with the flux 
hmited 10kspec+30kZADE sample of tracer galaxies. The aperture is deflned by the distance to the 5th nearest neighbour projected within 
zb 1000 km to redshift of the individual grid points. All three axis, distance from the survey centre in right ascension, distance from the 
survey centre in declination and distance (corresponding to redshift), are expressed in comoving h~-'^Mpc. The reconstructed overdensity 
field yields structures in a large range of comoving scales covering a spectrum of different overdensities at all redshifts reliably probed by 
the 10k zCOSMOS survey. The colour scale of the 1 -\- dp values is presented below each cone covering 500 h~^Mpc. The positions of 
galaxies are marked with yellow dots. 
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Fig. 12. — Distribution of the overdense structures in the 10k zCOSMOS overdensity field. The structures are defined as isosurfaces 
enclosing regions with l-\-6p > 1.5, 3, 5 and 10 (l + (5 > 1.8, 5.2, 11.2 and 31.6) from the left to the right, respectively. Comoving distances 
in the radial direction correspond to 0.1 < z < 1. The distribution of overdense structures matches the richness and complexity of the local 
cosmic web up to z ~ 1. The structures are delineated from the overdensity field reconstructed following the same scheme as in Figure [TTI 
All three axis, distance from the survey centre in right ascension (left hand side axis), distance from the survey centre in declination (right 
hand side axis) and distance (corresponding to some redshift; vertical axis) are expressed in comoving h~-'^Mpc. The figures are compressed 
by a factor of ~ 3.5 in redshift direction. 
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Fig. 13. — The 10k zCOSMOS underdense regions (voids). The structures are defined as isosurfaces enclosing regions with 1 -\- 6p < 
0.15 (6.67 times smaller than p(z)) and 0.25 (4 times smaller than p(z)) in the left and the right cone, respectively. Comoving distances in 
the radial direction correspond to 0.1 < z < 1. The structures are delineated from the overdensity field reconstructed following the same 
scheme as in Figure [TTI All three axis, distance from the survey centre in right ascension (left hand side axis), distance from the survey 
centre in declination (right hand side axis) and distance (corresponding to some redshift; vertical axis) are expressed in comoving h~-'^Mpc. 
The figures are compressed by a factor of ~ 3.5 in redshift direction. 
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Fig. 14. — Fractions of volume, galaxies, stellar masses and B-band luminosities enclosed in the structures defined by the given overdensity 
1 + (5p values, going from the left to the right. These fractions are presented in a narrow redshift slices = 0.2, starting from z = 0.2. 
The resulting fraction are dominated by cosmic variance. However, it is clear that the overdense structures are much more important in 
terms of their baryonic content than in the terms of the volume which they occupy. Note the different scale of the y-axis in the first panel. 




Fig. 15. — Comparison of the overdensities reconstructed at the positions of all zCOSMOS galaxies and those galaxies residing in the 
virialised structures. The histogram of the overdensities at the positions of zCOSMOS galaxies is presented in magenta. The shaded 
histograms represent the distribution of overdensities centred on galaxies which reside in the virialised structures: X-ray clusters (green 
shaded histogram), the optical groups with at least 2 detected members (black shaded histogram) and the optical groups with at least 4 
detected members (red shaded histogram). The histograms overlap each other. 
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Fig. 16. — Comparison of the zCOSMOS density field to the estimates of the bound structures. On the left and right figures is presented 
the zCOSMOS density field at 1 + (5p > 3 obtained with the 10k+30kZADE flux limited sample and the apertures deflned by the distance 
to the 10th nearest neighbour projected within ± 1000 km s~^. On top of the overdensity fleld are overlayed X-ray detected clusters (red 
circles, on the left) and optical groups with at least three detected members (magenta circles, on the right). Each circle is centred at the 
position of the deflned structures. Radii of the circles are scaled in RA — DEC plane to the X-ray luminosity for the X-ray clusters and to 
the number of detected objects in the optical groups. 
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Fig. 17. — The overdensity field of galaxies projected in redshift slices of 0.025 width, covering the redshift interval from 0.2 to 1. 
The overdensity field is reconstructed using the flux limited sample of galaxies within the apertures deflned by the 5th nearest neighbour 
projected within ± 1000 km s"-*^. Only structures deflned by the isosurface contours of 1 + (5p > 6.67 (red) and I -\- Sp < 0.15 (green) are 
shown. As a comparison, the optical groups in the same redshift bins with al least 3 detected members are overplotte d a s circles and X-rays 
clusters as squares. The sizes of the symbols to mark the positions of the virialised structures are scaled as in Figure [T6] The redshift slices 
0.3 — 0.375, 0.675 — 0.75 and 0.875 — 1 are dominated by the overdense structures. Large, RA — DEC-eyitended underdense structures are 
detected in the 0.45 - 0.5, 0.525 - 0.55, 0.575 - 0.6, 0.625 - 0.675 and 0.8 - 0.825 redshift slices. 
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Fig. 17. — Continued. 
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Fig. 17. — Continued. 
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Fig. 17. — Continued. 
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Fig. 18. — Comparison of the l + 6p =3 isosurfaces in the overdensity field reconstructed using the 10k zCOSMOS sample of flux limited 
tracer galaxies and the equivalent sample of galaxies in the mock catalogues. First left hand flgure: overdensity fleld reconstructed using the 
real 10k+30kZADE zCOSMOS galaxies. Rest of flgures: overdensity fleld reconstructed using the mock 10k+30kZADE samples. In every 
point the overdensity fleld is reconstructed by counting the number of objects with an integrated ZADE-modifled probability distribution 
in the apertures deflned by the distance to the 10th nearest neighbour projected within ± 1000 km s~-'^of the redshift of the grid point. The 
grid is regular in ARA = ADEC = 2 arcmin and Az = 0.002. The overdensity values are presented without the edge correction. The high 
complexity of the cosmic web is noticeable in all flgures. It appears that there are more structures on large scales (~ 10 h~-'^Mpc), visible 
above z ~ 0.8, in the real data than in the mock catalogues. To highlight the structures in the flgures, we omit the axis. The transversal 
axis are RA and DEC, covering the zCOSMOS area ~ 1 deg^, and the vertical axis is redshift in the range 0.2 < z < 1. 
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Fig. 18. — Continued. 
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Fig. 19. — Fraction of the volume with the overdensity 5p above a given value. The aperture in the overdensity field is defined using 
the distance to the 10th nearest neighbour projected to the redshift of the grid point within ± 1000 km s"-*^. A mock value is obtained by 
averaging results obtained from the individual 12 mock catalogues. The error is calculated as the standard deviations from the individual 
mock catalogue results. The continuous lines are for the real data, the dotted lines are for the corresponding mock catalogues. Volumes 
are estimated in [h~-'^Mpc]'^. Left: Statistics in 0.2 < z < 1. The dashed red lines are for the real data, the continuous blue lines are the 
corresponding mock catalogues. Right: Statistics in four redshift bins: 0.2 < z < 0.4 (green), 0.4 < z < 0.6 (magenta), 0.6 < z < 0.8 
(cyan) and 0.8 < z < 1 (black). 
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Fig. 20. — A comparison of the sizes of the structures above 1 + 5p > 3 in the data [red: V /Vmax smoothing to obtain N(z), magenta: 
smoothing equivalent to the smoothing applied to the mock N(z)] and in the individual mock catalogues (black curves). The curves 
correspond to the volume fractions contained within the structures of at least the size indicated on the x-axis, where size is measured in 
[h~-'^Mpc]^. There is not a single mock catalogue which contains as much volume as the real data in the large structures in 0.2 < z < 1. 
When we limit our statistics to 0.3 < z < 0.9, the data results falls within the statistics outlined by the mocks. Difference between the 
data and the mocks at the large sizes is dominated by the structure at 2; ~ 0.9. 



